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The competition between drift wave and interchange physics in general 
E-cross-B drift turbulence is studied with computations in three dimen- 
sional tokamak flux tube geometry. For a given set of background scales, 
the parameter space can be covered by the plasma beta and drift wave 
collisionality. At large enough plasma beta the turbulence breaks out 
into ideal ballooning modes and saturates only by depleting the free en- 
ergy in the background pressure gradient. At high collisionality it finds 
a more gradual transition to resistive ballooning. At moderate beta and 
collisionality it retains drift wave character, qualitatively identical to 
simple two dimensional slab models. The underlying cause is the nonlin- 
ear vorticity advection through which the self sustained drift wave tur- 
bulence supersedes the linear instabilities, scattering them apart before 
they can grow, imposing its own physical character on the dynamics. 
This vorticity advection catalyses the gradient drive, while saturation 
occurs solely through turbulent mixing of pressure disturbances. This 
situation persists in the whole of tokamak edge parameter space. Both 
simplified isothermal models and complete warm ion models are treated. 



PACS numbers: 



52.55. Fa — Tokamaks 
52.35.Ra — Plasma turbulence, 
52.65 — Plasma simulation. 
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I. Introduction — More Than One Eigenmode in Turbulence 



Edge turbulence in tokamak flux tube geometry has been treated by models working 
from a drift wave [1,2] or resistive ballooning [3,4,5] paradigm. Beyond computing edge 
transport, the main purpose of these computations is to understand the underlying physical 
character of the turbulence: drive and saturation mechanisms and free energy transfer 
channels. 

The pure situation of drift wave turbulence is slab geometry in a sheared magnetic 
field, and all of its basic processes can be found in the simplified two dimensional models in 
which it was first studied [6,7,8,9,10,11]. The idealised state is one of adiabatic electrons, 
towards which the electron dynamics parallel to the magnetic field B tends: 

n e eB • V0 — > B • Vp e (1) 

(where equality means the electrons are adiabatic, and where is the electrostatic potential 
and n e and p e are the electron density and pressure, and e is the electronic charge). This 
is the adiabatic response, and if it carries all the way to adiabatic electrons there can be 
no free energy release, because the average transport is zero if there is no phase shift in 
the direction of the electron drift given by — BxVp e . The pressure force B • Vp e in the 
Ohm's law conserves energy against the compression of parallel currents along field lines, 
V - ( J||B/S), which provides the back reaction in the electron pressure. Together, these two 
terms constitute the adiabatic coupling mechanism, which allows the free energy liberated 
from the background gradient to enter the parallel current. The more familiar Alfvenic 
coupling with B ■ V0 then allows the free energy to enter the ExB flow eddies. This is the 
free energy pathway between the background gradient and the turbulence which gives rise 
to drift waves, and collectively the adiabatic and Alfvenic couplings constitute the drift 
wave coupling mechanism p e <-> Jn <-> among the respective disturbances. It is important 
to note in this context that the parallel gradient cannot vanish for finite sized disturbances 
on closed magnetic flux surfaces (the flux tube geometry is globally consistent [12]), so this 
adiabatic response is always excited by electron pressure disturbances. 

The pure situation of resistive ballooning mode turbulence or ideal ballooning in- 
stabilities is an MHD (magnetohydrodynamics) model [13]. In this case, the adiabatic 
coupling mechanism is specifically neglected so that dissipative Alfvenic activity merely 
damps e<j)/T e towards zero instead of towards p e /p e (here, T e is the electron tempera- 
ture in energy units, and we assume p e = n e T e ). The only way that free energy can be 
transferred to ExB motion from the background pressure gradient in an MHD model is 
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through the magnetic curvature and gradient, which give rise to a finite compressibility of 
the diamagnetic current and hence the well known interchange effect. In toroidal geome- 
try, this leads to the ballooning mode since the interchange effect is destabilising only on 
the outer part of the torus, where the radius of curvature and gradient-B vectors become 
parallel to the pressure gradient. The ideal form of this dynamics as an instability in 
toroidal geometry was discovered in attempts to explain curious observations of magnetic 
fluctuations localised to the outboard side of the PLT tokamak [14]. These fluctuations 
were explained as ideal ballooning modes, which were treated conventionally [15] as well as 
with what became known as the ballooning transformation [16,17,18]. Resistive ballooning 
[19] was invoked to explain the existence of tokamak edge turbulence at pressure gradient 
values below the ideal stability limit and at amplitudes seemingly unattainable by other 
mechanisms [20,21], and has also been the focus of many linear MHD calculations [22]. 
More recently [23] , the ballooning paradigm was extended to situations including the two 
fluid Ohm's law (i.e., the adiabatic response), leading to the ballooning mode approach 
to edge turbulence [3,4,5], coinciding with the development of the earliest treatments of 
flux tube geometry, which were originally constructed with ballooning modes specifically 
in mind [24,25]. 

The new phenomenon which arose with the three dimensional flux tube models is the 
simultaneous presence of these two quite different types of eigenmode in the same dynami- 
cal system. Obviously, there is only one electrostatic potential, so it is not germane to treat 
both processes separately and simply add up the resulting growth rates and hence mixing 
length transport estimates (which one could do for linear eigenmodes if the wavenumbers 
and frequencies were significantly different) . We have a situation of ExB turbulence in the 
presence of a background pressure gradient, which is the basic free energy source for both 
eigenmode types. What differs is the energy transfer channel, that is, parallel dynamics 
or interchange forcing. The MHD model contains the interchange mechanism but neglects 
the drift wave mechanism by omitting the adiabatic coupling. The slab drift wave model 
neglects the interchange mechanism by omitting the gradient-B and magnetic curvature 
effects. One must obviously retain these interchange effects in a toroidal model, so the 
principal consideration becomes the adiabatic response: the neglect of that limits the va- 
lidity of the MHD model to situations in which the disturbances satisfy p e /p e <C e0/T e (cf. 
Eq. 1). For purely parallel (shear Alfven) dynamics this is satisfied if k±p s <C 1 [26], where 
p s is the drift scale at which simple linear drift waves principally occur. However, for ExB 
motion driven by a pressure gradient, the typical dynamical frequencies are comparable 
to cu*, the diamagnetic frequency, and hence ExB advection of the pressure will render 
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Pe/Pe ~ e0/T e , which is the situation of drift waves (note that this does not mean adia- 
batic electrons, as the phase remains arbitrary and the equality need not be very close). 
The MHD paradigm is therefore not useful in addressing this situation; moreover, it is 
impossible even to discuss drift waves within MHD due to the missing adiabatic coupling 
mechanism. 

The competition between the two types of eigenmode is strongly complicated by non- 
linear physics, when one addresses turbulence. For linear drift waves in the presence of 
interchange forcing, it can be shown that the interchange forcing can be the principal 
mechanism leading to a finite parallel current divergence and therefore departures from 
the adiabatic state, hence phase shift, over a wide range of parameters. This phenomenon 
is what gave rise to the "drift resistive ballooning" paradigm, in which it was originally 
attempted to order p s small in the dynamics [23]. Indeed, the scales of motion in that 
paradigm are derived from the assumption that interchange forcing balances the paral- 
lel current divergence. On the other hand, it was already known that while the linear 
polarisation drift by itself could only cause a weak nonadiabatic response, the nonlinear 
dynamics greatly transforms the mode structure of the turbulence [7], leading to robust 
turbulence even in models with no linear instability at all [10]. Early attempts to treat 
simple interchange forcing within such models showed that under certain circumstances 
the interchange forcing did not affect the turbulence even if it could cause significant lin- 
ear instabilities [27], a result which was confirmed by three dimensional studies [1]. If 
the nonlinear polarisation drift is strong enough to overcome the interchange forcing in 
the turbulence, it will provide the balance with the parallel current thereby maintaining 
electron nonadiabaticity, and the principal scale of motion will vary mostly with p s even 
though it is larger than p s . This is even true in two dimensional slab turbulence, where 
there is also dependence on dimensionless parameters such as the collisionality [10]. In 
three dimensional electromagnetic computations, a finite plasma beta was found to en- 
hance the nonlinear drift wave as well as the interchange (ideal ballooning) effects [2]. It 
therefore becomes crucial to resolve p s in any computation of tokamak edge turbulence, 
and crucial further to use a formulation of the magnetic geometry which represents both 
slab and toroidal mode structure equally well. For the question of what type of mode 
structure obtains is closely related to which of the polarisation or diamagnetic currents 
provide the balance with the parallel current under the total divergence. 

In this Article we treat the competition between drift wave and interchange physics 
for control of the turbulence in detail. Firstly working within the simplest model treating 
both eigenmode types (isothermal electrons, cold ions, where the difference to MHD is in 



4 



a single pair of terms) we choose parallel and perpendicular scales typical of edge profiles 
in modern tokamaks and sweep the plasma beta and collisionality parameters, which is 
the same thing as varying temperature and density at fixed gradient scale length. We 
compare the details of the mode structure — amplitude ratios, parallel structure, phase 
shifts, cross coherence, and drive spectrum — to control cases reflecting the pure situation 
of both eigenmode types. We also directly measure the details of the vorticity dynamics 
- mean squared divergences of the pieces of the total current and mean squared transfer 
effects — to find which two mechanisms are in statistical balance. This gives a rather more 
complete picture than one would have merely with the usual time traces of transport levels 
and contour plots showing spatial morphology. 

We find that the transitions to both ideal and resistive ballooning turbulence are well 
defined, and that the "working area" of tokamak edge density/temperature parameter 
space is occupied by drift wave turbulence. By studying situations where an interchange 
driven linear eigenmode makes the transition to fully developed turbulence, we also show 
that the reason for the drift wave character is the nonlinear polarisation drift, which 
replaces the diamagnetic current as the main balance for parallel currents and nonadiabatic 
electron activity. It is always dangerous to form one's insight for turbulence from the 
properties of linear eigenmodes, and this error is essentially what happened to the drift 
resistive ballooning paradigm. The conclusions of linear theory are not valid for turbulence, 
because drift wave turbulence is very much more robust than linear drift waves are: in a 
variant of the well known ExB shear suppression scenario [28] , the rms vorticity level of the 
turbulence is larger than the ideal interchange growth rate in this steep gradient regime, 
and so the small scale linear instabilities arising from ballooning effects are superseded 
by the nonlinear drift wave instability which imposes its own physical character on the 
dynamics. Finally, though the subject deserves its own treatment in a separate paper, 
we address the question of the saturation mechanism, with the result that the vorticity 
nonlinearity always drives while the pressure one provides saturation, as in the pure drift 
wave case. 

II. The DALF3 Model 

The most basic model of ExB drift turbulence containing both drift wave and inter- 
change physics is a four field model in tokamak flux tube geometry which we call DALF3, 
for three dimensional drift Alfven turbulence [2]. The magnetic geometry is assumed to 
be a set of nested, toroidal, axisymmetric flux surfaces, on which one can define a glob- 
ally consistent flux tube [12]; this is a matter of defining Hamada coordinates aligned to 
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the equilibrium magnetic field B for which there is only one nonvanishing contravariant 
component, which is a flux function, and applying the parallel boundary condition after 
one connection length Ly, ensuring that all degrees of freedom in the flux tube actually 
exist on the flux surface for which it is a model. A technique called "shifted metric" then 
re-adapts the geometry so as to represent slab and toroidal mode structure equally well, 
by providing for a perpendicular metric which is locally diagonal within each drift plane 
perpendicular to B [29]. Following the local drift scale ordering, the dependence of the 
metric on all coordinates except the one parallel to B is neglected (cf. Ref. [25]). 

The dynamics consists of the fluctuating electrostatic potential (0) and electron pres- 
sure (p e ) as state variables, those acting as force potentials, and the parallel current (Jy) 
and parallel ion velocity (uu) as their corresponding flux variables. Their set of equations 
is formed by the two conservation laws for charge (V • J = 0) and electron thermal energy, 
and then the parallel component of the equations of motion for the electrons and ions. The 
parallel magnetic potential (A\\) is given by Jy through Ampere's law, and then cf> and A\\ 
serve as stream functions for the ExB velocity (ve) and the magnetic field disturbances 
(Bj_), respectively. Traditional flute mode ordering is assumed, with parallel and perpen- 
dicular components satisfying L^ 1 ~ Vy <C L^ 1 ^,Vj_, where L± is the perpendicular 
scale of the background. The background is a set of constant parameters except for where 
ve ■ V or Bj_ • V act on the background pressure gradient. The only nonlinearities are 
ve ■ V and Bj_ ■ V, the latter which is part of Vy. The derivation is standard; one can 
consult Ref. [2] for an example (this model is derived from that one by setting T e to zero). 
With the neglect of the temperatures the model is not different from the one by Braginskii 
[30] under the local drift approximation [31], retaining electromagnetic parallel dynamics 
[32]. 

A more general model will be briefly invoked in the Section IX to show that the results 
are not dependent on the neglect of the temperatures, but the usefulness of DALF3 is its 
closeness to resistive, reduced MHD: the two models differ only in the presence or absence 
of the adiabatic coupling mechanism and of the (nearly negligible) diamagnetic heat flux 
in the pressure equation. This allows investigation of the special effects of the adiabatic 
p e <-> J|| coupling in direct fashion. 

The equations are normalised in terms of drift scales p s and c s /L±, defined by 

2 _ C 2 MjT e 2 _ Te /r>N 

Ps ~ ~~e 2 B 2 ~ Cs ~W z 1 J 

This is the standard gyro-Bohm normalisation, leading to a diffusion level Dgb = p 2 s c s /L±_ 
and a flux velocity level Tqb = c s (p s /L±) 2 . The dependent variables are scaled in terms 
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of ecf)/T ei p e /p e , J\\/n e ec s , and u\\/c s , and the auxiliary one in terms of A\\/j3Bp s , where (3 
is one of the parameters discussed below. The dependent variables are further scaled with 
an additional factor of p s /L±, such that a normalised (j) = 1 refers to an unnormalised 
e(j)/T e = p s /L±. The parallel gradient and divergence are normalised against L||/27r, while 
the perpendicular gradient is normalised against L±. The normalised equations appear as 

^lvi?=BV|,i-;C(p e ) (3) 
~cL4.ii dJ\\ _ ~ ~ 

j t (p e + p e ) = SV|| " /C(p e - 4>) (5) 

e -^- = - V ||(Pe+Pe)+//||V||U|| (6) 



with the Ampere's law 



Jn = -V\A* (7) 



serving as a constitutive relation (since the eigenvalues of — are positive, Jy and A\\ 
are in essence proportional and can be used interchangeably as dependent variables). Note 
that in the DALF3 model p e = n e (we thereby avoid the misconception that drift wave 
dynamics is driven solely by the density gradient — cf. Ref. [10]). 

For the present purpose we use the simplified form of the sheared flux tube geometry, 
neglecting all finite aspect ratio effects except for the existence of the interchange operator 
/C. The factors of the magnetic field strength in Eqs. (3,5) are set to unity. The connec- 
tion length L|| becomes 2ixqR, with q the standard field line pitch parameter and R the 
major radius. We use the shifted metric coordinate system [29], which uses a different 
y-coordinate at each location in s, labelled yu for the particular field aligned coordinate 
system which has a locally orthogonal perpendicular Laplacian operator at s = Sf~. When 
taking derivatives in s we then have a shift in the y-coordinate, since the various yk are 
different. The perpendicular differential operators for the k-th coordinate system at s = 
are 

d 2 8 2 „ ( . d d 



V ± = iro + 7T2 £ = U B [ sins-- +coss-^ ) (8) 

ox 2 ayf. \ ox oy^ J 

where wg is a parameter, and the centered partial s-derivative at s = Sk is given by 

df , s f(x,yk+i - Q!fc + i,Sfc + i) - f{x,yk-i - afc-i,Sfc-i) / n >. 

— {x, y k - a k ,s k ) = (9) 

as s k +i - Sk-i 



7 



with shifts given by 

a k = ss k x (10) 

and s parameterising the magnetic shear. We note that in each coordinate system B- Wy k = 
for the equilibrium B, because B- Vx = 0. We note as well that the shifts do not affect the 
partial y-derivative, which is hereafter written without the subscript k. Similar centered 
differences are defined at the mid-node positions s k+ i/2, which shift y k and y k +i onto 
Vk+i/2, as required by the second order upwind numerical scheme (cf. Ref. [33], which 
adapts the second order MUSCL scheme [34,35] to this problem). Accordingly, the entire 
set of geometric quantities (metric, magnetic field strength and unit vector, interchange 
operator) are defined at both nodes and mid-nodes k and k + 1/2. It is important to 
note that the shifting does not merely reverse the field aligning, as with each s k a global 
constant, we still have B y — since each a k is strictly a function only of x. 

The advective and parallel derivatives are given by 

| = | + v £ .V V„=b.V (11) 

where ve and b have contravariant components given by 

•*-£ ^4 < i2) 

The appearance of the background p e under derivatives represents the gradient drive terms, 
through 

p e = -UJ p X (13) 

where uj p is a parameter (note choosing u p = 1 is equivalent to choosing L±_ = L p ). 

The dependent variables are defined on a domain 

n n n n 

< x < — — - — <y<— -TV<S<TV (14) 



AK AK K K 

where K and A are parameters. The boundary conditions are Dirichlet for x (dependent 
variables vanish at the boundaries), periodic for y, and in s they satisfy 

f(x,y k+N ,s k+N ) = f(x,y k ,s k ) (15) 

where iV is the number of nodes, noting that the two points are on differing coordinate 
systems, as in the parallel derivatives. The profiles are maintained by boundary source 
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terms, damping the k y = component of p e by a Gaussian of amplitude 0.1 and a 1/e 
width of 0.1 times the domain width, centered at each boundary. 

The physical parameters are 

2 2 2 

h 47rp e / qR\ ^ m e f qR\ L ± m e f qR\ 

for the dissipative kinetic shear Alfven dynamics, and 

L± 2L ± . (qR\ 2 

for the gradient drive, interchange forcing, and the sound waves, with an arbitrary 
parallel diffusion coefficient. Here, T ue is the standard Braginskii collision time for the 
electrons [30] . Nominal parameters corresponding to a typical tokamak plasma edge are 

(3 = 1 fa = 5 C = 7.65 u B = 0.05 e = 18350 s = 1 (18) 

with u p = 1, roughly reflecting physical parameters: 

n e = 3 x 10 13 cm" 3 T e = 70 eV B = 2.5 T 

(19) 

#=165 cm L ± = 4.25cm q = 3.5 

The grid for all cases consists of 64 x 256 x 16 equidistant nodes in {x,y,s}, with a 
domain size of 207r x 80-7T x 2-7T, respectively (K = 0.025 and A = 4). The time step is 
0.05, and the initial state is a random phase distribution for p e , with amplitude spectrum 
[1 + (^/vT28ii' 2 ) 4 ]" 1 normalised to an rms amplitude of 3.0, and with <f> = A\\ = u\\ = 0. 

III. DALF3 Energetics 

The physics of turbulence cannot be assessed from transport scalings or pictorial 
morphology; quantitative diagnostics of the amplitude, energy, and energetics spectra, as 
well as several cross-variable statistical measurements relating the transported quantities 
to the ExB flows that transport them, are needed [1,2,10,11,27,33,38,41]. In the present 
case, the diagnosis of the energy transfer channels between the equations for p e and (p are 
among the more important. 

Since the energetics of the DALF3 model will form a central part of the analysis, 
we review it once more, within the treatment of magnetic coordinates just outlined. To 
form the energy theorem we multiply Eqs. (3-6) by —(f), Jy, p e , and u\\, respectively, and 
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integrate over the spatial domain, assuming that total divergences vanish. We find the 
following, in which the integration operation is denoted by the angle brackets, 



d_ 
di 



1 1 

2S2 



V 



JhVi 



p e IC(f) 



§t (l J \\ (Ml + ) = V||p e ) - ( J|| V|,0) - C 



Jll 



dt \2 lPel 
d_ 

di 



- \ J \\^\\Pe) + («|| V||P e ) + (Pefccp) + U p ( -p, 



0(P 

dy 



2 hi 



- («||V||p e ) ( U\ 



My 

ay 



^ii 



Viiuii 



(20) 
(21) 
(22) 

(23) 



The quantities under the partial time derivatives are the ion drift energy, magnetic and 
electron kinetic energy, thermal free energy, and sound wave energy. Terms appearing in 
the right sides of two equations with opposite sign represent the transfer channels, while 
isolated terms reflect the sources and sinks. The only important source is the one in the 
pressure equation, ExB advection down the gradient. Currents flowing down the gradient 
cancel out of the energetics, and sound wave motion down the gradient is always negligible. 
The sinks are resistive friction (C), sound wave viscosity (/iy), and the subgrid dissipation 
resulting from the nonlinearities, especially v E ■ Vp e , transferring free energy to arbitrarily 
small scales — out of the spectrum on the high-/c^ side. 1 

There are two ways for free energy to be channelled from the source in the pressure 
equation into the ExB turbulence: the drift wave mechanism, which couples p e <-> Jy «-> 
through the parallel dynamics (Vy), and the MHD compression mechanism, which couples 
p e <-> <p directly through the interchange forcing (/C). One of our tasks will be to determine 
which of these is the more prominent. If we remove the /C terms only the drift wave 
mechanism is present; if we remove the Vy (p e +p e ) term in the Ohm's law (Eq. 4) and the 
VyJy term in the pressure equation (Eq. 5) only the interchange mechanism is present. 
These form the two relevant control cases. The linear transfer effects (Vy and JC) would 
proceed unidirectionally for linear or coherent nonlinear systems, but with the vg- V terms 
providing a quasirandom character the transfer effects should be understood as exchange 



1 Conventional treatments treat this through artificial viscosity terms; herein, the sub- 
grid dissipation is an integral part of the numerical scheme; in any case, it is a defining 
feature of high Reynolds number turbulence, with the total dissipation independent of the 
size or form of the dissipation operator. 
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mechanisms — as we will see, their rms levels are often much greater than their mean 
values. 



Half squared amplitude, free energy, and source and sink profiles and spectra can be 
defined, for example, 

1 v 1 - 1 i ~ 2 
M k v) = -^Ylj^.Yl 2 \<t>( x i> k y> s k) (24) 

for the half squared amplitude of <p, with k y = IK for the Z-th Fourier component, and the 
iV's giving the number of grid nodes in each dimension. Half squared means simply 1/2 of 
mean squared; the factor of 1/2 is used because some of the half squared amplitudes (e.g., 
A p for p e ) also appear as free energy components. Further details of how these quantities 
are defined and measured can be found in Refs. [1,10]. 

IV. Drift Wave and Interchange Control Cases 

Here we re-examine two dimensional drift wave and interchange turbulence models, 
in order to provide a control for the analysis below. From the DALF3 model we neglect 
the s-direction, and set /3 = p, = 0, substituting the resulting collisional Ohm's law into 
Eqs. (3,5). We replace — C _1 V| with the positive constant parameter D for the dynamics 
involving Jy, and neglect sound wave motion. This reverts to the well known Hasegawa- 
Wakatani model [7] with interchange forcing. We also add a damping term 7^ to allow 
the interchange case to saturate. The resulting two dimensional control model is 

±V 2 J = D{4> -p e )- UB ^- 7d (25) 

-^(Pe+Pe) = D(4>-p e ) -^V~Q^ - W B-^(Pe-(p) (26) 

The drift wave version uses D = 0.05 and ub = 7d = 0, and the interchange version uses 
D = 0, ujb = 0.05, and 7^ = 0.01. For the drift wave model we eliminate the variations of 
the k y = components at each time step. Within the single drift plane, we use the same 
domain size and initial state as in the three dimensional model. We start both models 
with the given initial state and run to t = 1000. The time evolution of the half squared 
amplitudes of (f> and p e are shown for both cases in Fig. 1, demonstrating saturation. 

These two models differ fundamentally in the degree of coupling between p e and 0, 
and we will now see that this produces an unequivocal distinction in their mode structure. 
For the saturated stage 500 < t < 1000 we show the average cross coherence and the 
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phase shift probability distribution spectrum between p e and <f> for both models. For the 
cross coherence we sample p e and at each node in the domain and at 50 equidistant 
time intervals, normalise each to its standard deviation, and plot the contours of the 
resulting histogram. For the phase shifts, we sample a p< j,(l) = lmlog(p ei *0/) for each 
Fourier component k y = IK, at each node in x and the same 50 equidistant time intervals, 
and plot the contours of the spectrum of resulting histograms. The signature of the drift 
wave coupling mechanism is very clear: For the drift wave case the cross coherence is 
significant and the phase shift distributions are narrow and close to zero (essentially all 
activity below about a v $ = 7r/4). It is important to note that this affects every mode in 
the spectrum, not just the dominant drive range of k y p s ~ 0.2. For the interchange case 
we find the opposite situation of nearly no cross coherence and phase shifts close to n/2 
since every mode is strongly driven. 

This is the mode structure information against which we will compare the results of 
the three dimensional computations. On the other hand, the drive spectra are very similar 
and it is important to note that these and other information relying principally or solely 
on the properties of n e or p e alone cannot distinguish these two eigenmode types. 



We now run computations within the full DALF3 model and compare them to the 
resistive MHD model at the same parameters. The MHD model is a reduced MHD case, 
in which under the drift approximation the ExB flow is dynamically incompressible, and 
so the perpendicular electric field is electrostatic, just as in the drift wave model. The 
pressure is the total pressure in MHD, but within a cold ion model the total pressure is 
simply p e . The cold ion model without temperature effects is therefore the closest to the 
corresponding MHD case. The only difference in the physics is the adiabatic coupling 
mechanism p e <-> Jy, represented by the terms V||p e m the Ohm's law and Vy J\\ in the 
equation for p e . In the MHD model this coupling is absent, and the diamagnetic flux 
divergence is neglected in the pressure equation; in all other respects the DALF3 and 
MHD models are identical (in MHD the distinction between warm or cold ions is lost 
through the neglect of the diamagnetic flow vorticity and gyroviscosity, and the neglect of 
the adiabatic coupling mechanism in the electrons). The MHD model appears as 



V. Direct Comparison of DALF3 to Resistive MHD 



J_d_ 

B^dt 



Vi0 = SV|| 



(27) 




V||0-CJ|| 



(28) 
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with the Ampere's law, 

J\\ = -Vil|| (31) 

All other considerations are as in Section III. Simple comparison to Eqs. (3-6) shows that 
the resistive MHD model is an entirely included subset of the DALF3 model. 

We run both models with the fixed parameter set shown above, and vary C from very 
low to very high values. The cases shown are for v = {1,2,5,10,20,50,100,200}, with 
C = 2.55 x v, and jtiy = (allowing the numerical scheme to provide dissipation at the 
very shortest parallel as well as perpendicular wavelengths). The normalised v is L±/c s r ue , 
and C = 0.51/iz/. 

The first quantities to show are the practical ones, the transport and the relative 
amplitudes between cf> and p e , which appear in Fig. 3, with the four main drift wave 
cases with v = 1, 2, 5, 10 checked with doubled resolution in the drift plane. If the MHD 
model predicted these correctly the DALF3 model would represent an unnecessary level 
of complication. We find that the transport scaling does in fact disagree, strongly. We 
also find that the rms amplitude ratios <fi/p e disagree. In drift wave dynamics, p e and cj) 
are held together by the dissipative adiabatic response, but in MHD this Alfven dynamics 
dissipates <j) only. Since MHD neglects the electron pressure in the Ohm's law, it implicitly 
assumes that <j> 3> p e in any situation in which the geometry sets the parallel scale as 
it does here. But with only damped by the Alfven effects, the MHD model actually 
delivers a result contrary to its validity, that p e > <fi. The DALF3 model finds the result of 
the adiabatic response, namely, 4> ~ p e . Both models give a power law for this amplitude 
ratio at low C. The position at which these two lines cross defines the regime boundary 
between the two models. For C larger than this critical value, both models find <j> > p e . 
It follows that the resistive MHD model is valid only for C beyond this value. Because 
the interchange forcing enters through the parameter ub, we find the boundary for the 
transition to resistive ballooning turbulence is at 

u B = Cu B = 1.02-£- (^) ' > 1 (32) 

-tiT ue \V e J 

The transport curves both reach maxima and would eventually converge for very large C; 
this is due to the fact that both approach two dimensional interchange dynamics moderated 
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only by sound waves and nonlinear \ E ■ Vp e cascades for values of C so large that Jy is 
essentially negligible. Such values are however well out of the parameter range of interest 
for tokamak edge turbulence. 

Beyond the simple scalings, we turn to the mode structure in the transition range 
given by v = {5, 10, 20}. The most important diagnostics are the ones which tell us about 
the relationship between <p and p e , since this is where the models are different. In Fig. 4 
we examine the cross coherence between p e and <p f° r both models. The control cases 
appear in Section IV for drift wave and interchange turbulence. We find that the MHD 
model shows uncorrelated behaviour for all three values of v, while the DALF3 model 
shows moderate correlation even for the largest v. The adiabatic response coupling <fi to 
p e through J|| is always significant in this parameter range. 

Turning to the phase shift distributions for each k y component in the spectrum, we 
find the MHD model predicting the ideally unstable values of = n/2 for all k y , while 
the DALF3 model shows the transition between and tt/2 at low k y as v is increased. We 
see why we need these two diagnostics together: the cross coherence reflects the turbulence 
as a whole, while the phase shift spectrum shows each component. We can see that the 
transition is in the longer wavelengths, while the cross coherence is influenced also by the 
smaller scale activity which remains in a drift wave mode structure. In the drift wave 
regime, however, the entire spectrum functions as a unit, and to the extent the longer 
wavelengths cannot separate and form the transition to interchange dominated turbulence, 
the longer wavelengths are controlled by the shorter ones. 

For the fundamental reasons for this state of affairs we turn to the energetics. Recall 
that the only transfer mechanism between p e and <p in the MHD model is the interchange 
forcing. For drift wave turbulence in slab geometry this coupling occurs through the 
adiabatic response, mediated by Jy , and so the free energy liberated from the background 
gradient can only reach the ExB eddies through the parallel dynamics. In the DALF3 
model both coupling mechanisms are available, and so the question arises as to which is 
most important energetically as well as in the correlations. In Fig. 6, we measure the 
energy transfer terms affecting the ion drift energy Eq. (20), including the nonlinearity 
which operates between various wavelengths but within the ion drift energy. The relevant 
quantities are the rms levels of these transfer channels, since their importance to the 
turbulence is in their maintenance of each part of the energetics at finite levels, whatever 
the instantaneous sign. The results may be compared to those for drift wave and for 
interchange turbulence. 
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For the drift wave cases the nonlinear polarisation drift (v E -VV^) is largest, followed 
by the nonadiabaticity (Vy J\\), while the interchange forcing (ICp e ) is subdominant. It is 
especially important to note that the interchange forcing cannot account for the rms level 
of the nonadiabaticity, while the nonlinear polarisation drift is more than large enough. 
The nonlinear (v E ■ VVj_0) and linear (V\[d(j)/dt\) parts of the polarisation drift largely 
balance each other, with the excess maintaining the nonadiabaticity. 

For pure (MHD) interchange turbulence, especially at low collisionality where the 
resistive dissipation keeps <p small, the interchange forcing and parallel current divergence 
are in rough balance, and when <p and hence the polarisation drift is larger, the interchange 
forcing becomes larger than the nonadiabaticity. In the DALF3 model containing both 
eigenmode types, the transition from drift wave to interchange transfer dynamics occurs 
in the v > 10 range, at which all three effects are of comparable importance. 

We note now as in the end of Section IV the important fact that if we look only at 
the morphology of one dependent variable, p e for example, that we cannot easily see the 
difference between the two models. We can see this in the spectra and parallel structure, 
although if we have this data for both p e and <p the conclusion is obvious. We show the 
amplitude spectra for the same six cases, in Fig. 7. The spectra of p e look quite similar 
for the three cases v = {5, 10, 20} in both models, but the relationship to (p is different. In 
Fig. 8 the drive spectra are shown, and these too are very similar. Most importantly, the 
scale of motion is not a distinguishing factor between the two models. This means that 
experiments that measure the statistics or morphology only of density fluctuations, which 
are the most accessible in a hot plasma, can do little to ascertain what the dynamics of the 
turbulence actually is. 

We find a similar state of affairs in the parallel structure. The amplitude envelopes 
of the state variables are shown in Fig. 9, and the variation of the transport flux with s 
is shown in Fig. 10, each for the same six cases (for these parallel structure figures the 
k y = mode is stripped so as to concentrate on the part of the disturbances which lead 
to net transport down the gradient, noting that v% = —d(p/dy). Here we find that the 
outboard to inboard asymmetry of the density disturbance activity is also not a distin- 
guishing factor, as the mere presence of interchange forcing, even at the low levels implied 
by uj b = 0.05, is sufficient to make a difference in the overall dynamics. But this difference 
is a quantitative difference. The important part of the parallel structure is the much lower 
degree of asymmetry in h e compared to either <fi or p e in the DALF3 model. Although 
p e is strongly ballooned, the quantity which actually liberates free energy is actually /i e , 
whose relative lack of ballooning in DALF3 for ub ~ 1 reflects the importance of field line 
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connection and drift wave coupling. The fact that p e (or </>) is ballooned is actually not 
very relevant to the question of what type of turbulence this is. 

Taking all these diagnostics together, the qualitative difference of the DALF3 model 
with vb < 1 to the two dimensional slab drift wave model is negligible. Although the pres- 
ence of the interchange effects is noticeable, the omission of the adiabatic p e <-> Jy coupling 
leads to catastrophic changes. The basic mode structure features which are able to tell the 
difference between drift wave and interchange turbulence decide the contrast unequivocally 
for drift wave mode structure, and therefore dynamics, for the regime described by vb < 1- 

VI. The Role of Interchange Forcing in Drift Wave Turbulence 

It is instructive to examine the effect of the interchange terms, which do the in- 
terchange forcing, on simple linear drift waves with prescribed, constant k± and fey. This 
version of the linearised equations takes a single Fourier component by setting — > — k\, 
and Vy — > ik\\, and d/dy — > ik y , and for the interchange forcing K, — > iusky. To find the 
fate of the instabilities, all we need to know is the response of h e to cj> in order to find the 
dispersion relation, because the resulting phase shift is closely related to the linear growth 
rate. The vorticity equation and Ohm's law, Eqs. (3) and (4), are all we need. Using the 
relations cu* = uj p k y and u c = ivsky, their linearised form is 

ujk]_4>Vi = k\\ J\\ k - uj c p e k (33) 

w*) A|| k + (fkj + iC) Jj| k = -k\\ (p ek - <^ k ) (34) 

Eliminating J|| k , we find the equation for the response. It is simplest to do this in the 
collisional regime, where C is the dominant effect in Eq. (34) . We then find 

fcf - iCk\u „ 
^ = jfcj + lC u c ^ (35) 

As with drift waves proper, the consideration of the phase shift is only important when 
fey is finite, enough so that fe| is the largest term in both numerator and denominator of 
Eq. (35). If this is so, then the phase shift is small. It then follows that the eigenmode has 
the same properties as a drift wave, merely with a different excitation mechanism for the 
phase shift. But the role of the interchange forcing in this case is simply to excite a phase 
shift and thereby cause the turbulence to be driven by the background gradient. This is the 
same role the polarisation drift has for drift waves, and the same role its nonlinear version, 
the vorticity nonlinearity, has for drift wave turbulence. Even in this regime dominated 
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by parallel dynamics, the interchange effects can assume control and change the mode 
structure if they dominantly cause the phase shift. We note here that the fmiteness of 
fey for finite k y is guaranteed by the combination of magnetic shear and toroidal topology 
(closed flux surfaces), that is, field line connection [12]. 

VII. Linear versus Nonlinear Mode Structure 

The reason that the result of this investigation for the turbulence is so different than 
what one might expect by looking at linear instabilities is that the nonadiabatic character 
of drift waves, especially, is greatly changed by their associated nonlinearities. The one with 
the strongest effect is the nonlinear polarisation drift — equivalently, the self advection of 
unsteady ExB vortices, or eddies — which is the mechanism through which the turbulence 
has such strengthened parallel current dynamics, and hence more net nonadiabaticity of 
the electrons than those linear eigenmodes. This is the reason for and importance of the 
distinction between drift waves and drift wave turbulence. 

We can investigate this directly by starting the dynamics with an initial state in the 
linear stage with a linear instability strongly driven by the interchange forcing, and then 
tracking the evolution of the mode structure through the saturation stage, which is the 
transition from the linear stage towards the turbulence, and then into the fully developed 
turbulence in which there is no detailed memory of whatever transpires several correlation 
times before. 

The initial state is the same as before, just with an amplitude of ao = 3 x 10~ 10 
instead of the usual 3.0. We take a single case, that given by the nominal parameters in 
Eq. (18), but with v = 10 (hence C = 25.5 and vb = 1.25) and n\\ = 0. The resolution 
in the ^-direction is increased fourfold to diminish the effects of linear grid modes (the 
same domain size as before is used in all three dimensions, but the node count is now 
256 x 256 x 16; grid modes result from the failure to resolve the rational surface density in 
x for the largest k y component on the grid). The time traces of the amplitudes, transport, 
and energetics are shown in Fig. 11. The linear instability develops, and reaches a constant 
growth rate which is quite large, about 0.15 in the units of c s /L±_ (in which the ideal 
interchange growth rate is about ujb 1 ^ 2 = 0.22). When the disturbances reach finite 
amplitude, the growth rate rapidly drops to zero. This is the onset of saturation (the 
exact moment for this case is t = 201). Then, the growth rate begins to fluctuate, rises 
away from zero for a short while, and then approaches zero again up to some statistical 
fluctuation on the scale of a correlation time. This is the mode structure adjustment 
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stage, beginning for this case with the next rise of T T above zero at t = 212 and ending 
at t = 345 with the next drop of T T back to zero. The transport overshoots before 
settling, peaking near t = 200, quickly dropping by half to about t = 210, and then slowly 
strengthening throughout the adjustment stage until statistical saturation is reached after 
about t = 400 for the turbulence and after t = 600 for the larger scale ExB flows. For 
lower C the linear instability is confined to rather smaller scales, with more immediately 
efficient nonlinear decorrelation, and the overshoot is barely noticeable. In either case, the 
transport and overall growth rate reach statistical saturation only some several At = 100 
after initial saturation, by which time the drive and amplitude spectra have their final 
form. (Preliminary studies of these effects with only 32 x 128 grid nodes in the drift plane 
showed much stronger overshoot, with a longer period of statistical fluctuation before 
saturation — by comparison to the result shown here with full resolution one can conclude 
such behaviour is a manifestation of too few degrees of freedom in the turbulence.) This 
is the stage of fully developed turbulence, which we can see has a very distinct separation 
in time from the prior two stages. Full saturation also occurs only after the zonal flows 
(the k y = and k\\ = component of </>), to which the turbulence is energetically coupled 
[36,37], reach statistical equilibrium, in this case after about t = 600. 

These time traces show evidence for complete supersession of the linear instability 
by the drift wave turbulence. The last frame (lower right) of Fig. 11 compares the rms 
vorticity (mostly at 0.3 < k±p s < 1.0; see below) to the linear growth and nonlinear drive 
rates. In the conventional mixing picture of turbulence driven by a linear instability, the 
amplitude grows until the nonlinear rms vorticity O r ms (where O = |Vxv£;| = Vj_0) is 
comparable to the linear growth rate 7^. The eddies comprising the turbulence are then 
said to be driven by the instability and decorrelated nonlinearly by their vorticity, leading 
to the estimates of linear mixing length models: 

1l ~ k_LV k± -Dmix ~ ^ ~ -p 1 (36) 

where -D m i x is the turbulent diffusivity. In this situation, however, we have found not only 
that O rms is much larger than 7^, but also that the overall drive rate r + is significantly 
lower for fully saturated turbulence than either the drive or overall growth rates in the 
linear stage. This means that the turbulence is producing its own vorticity through the 
process of nonlinear self sustainment [10], and the reason that the linear instability at small 
scale has so little effect on the turbulence mode structure, as we have seen, is that small 
scale structures are scattered apart by the vorticity before they can grow. The dynamics 
therefore does not feel the linear instability in these scales. This is a generalisation of the 
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ideas of linear mixing length models as well as of the linear ExB shear-suppression scenario 
[28], simply instead of the instability being suppressed by background ExB vorticity (Qe = 
dYs/dx), it is suppressed by turbulent ExB vorticity: 



where -D S catter is the scattering diffusivity of the turbulence (nonlinear dissipation effect 
goes mostly through • Vp e ). This is the ultimate reason for the mode structure changes 
from the linear to the turbulent stages, which we now document. 

Fig. 12 shows the spatial morphology in the linear, saturation, and fully developed 
turbulent stages (only 1/4 of the y-domain is shown). Up until initial saturation, the 
properties of radial (Vx) interchange flows are pronounced, with the spatial morphology 
looking very much like that of buoyancy driven flow plumes (cf. t = 191). Initial saturation 
occurs when the nonlinearity is strong enough to radially displace the density structures 
and to break up the plumes (t £ 200) . The turbulence begins as both these structures 
and the ExB flows isotropise as a result of both main (v E ■ V on the state variables) 
nonlinearities (215^,^^,230). It then becomes established as the structures and eddies 
(coupled by the adiabatic response) expand in scale until the various linear and nonlinear 
transfer processes are in balance (400 ~ t ~ 600) . The isotropisation process is the same as 
that observed in the earliest collisional drift wave turbulence computations [7]. On top of 
the turbulence, the zonal flows develop as eddies serendipitously line up in the y-direction 
and then reinforce themselves by tilting the eddies into an energy-losing relationship vis- 
a-vis the flow [38]; these are visible as vertical stripes in the pattern of (j)(x,y). These 
striking visual changes provide motivation for the statement that the properties of linear 
instabilities have no a priori relevance to the turbulence. But to really find out whether 
the linear stage is relevant, we have to examine the mode structure in both stages. 

The diagnostics we have been using for the turbulence show important differences be- 
tween the linear, saturation, and fully developed turbulent stages. The dynamical transfer 
spectra for the ExB vorticity are shown for the transition stage in Fig. 14, noting that the 
state quickly approaches what we see in Fig. 6. In the linear stage, the nonlinear polarisa- 
tion drift is obviously negligible and there is a tight balance between interchange forcing 
and nonadiabaticity (which also shows that the linear polarisation drift, which involves 
the partial time derivative, is subdominant). As in the interchange turbulence cases, the 
linear stage of the DALF3 model shows a close balance between interchange forcing and 
the parallel current divergence. We therefore have an interchange-dominated linear insta- 
bility. Through the saturation stage, the polarisation drift enters, and in the end state the 
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polarisation drift has become the dominant cause of the finite Vy J\\ and hence electron 
nonadiabaticity h e = p e — 4> and hence the phase shifts necessary for finite free energy 
access and, ultimately, net transport, both of which are proportional to — § dy p e (d(f) / dy) . 
The parallel structure, displayed in Fig. 15, also shows interesting changes between the 
linear and turbulent stages. The degree of ballooning, especially in h e is markedly reduced 
in turbulence. Development of the nonlinear mode structure takes longer for this particu- 
larly three dimensional signal, following the slowness of the transport of the adiabatic part 
of the system along field lines, through sound wave dynamics. 

Finally, the phase shift distributions are examined in the stages of linear growth, 
nonlinear structural adjustment, and fully developed turbulence in Fig. 16. For this choice 
of parameters the linear instability is very strong (7^ m 0.15), a significant fraction of 
the nominal ideal interchange growth rate (7/ = ujb 1 ^ 2 = 0.22). Its dominant range is 
about 0.2 < kyp s < 1.0. When this range becomes nonlinear, the first consequence for the 
phases is that the large phase shifts in the range k y p s > 0.3 are eliminated, replaced by a 
wide distribution which also includes negative values. At late times the long wavelength 
component finds its structure — the frame on the right of Fig. 16 shows the same character 
as the corresponding case in Fig. 5. Indeed, the saturated state is independent of the 
previous history beyond a few correlation times; here, the correlation time is about 6, and 
long term, large scale memory persists no longer than about ten of these. With the short 
correlation time of the turbulence, the initial state with which the turbulence is reached 
is not relevant to consideration of its mode structure or underlying dynamics. It follows 
that this applies also to linear instabilities. Whether these are relevant depends on the 
particular case and can be demonstrated only by diagnostics taken on the fully developed 
turbulence in its own context. 

In the present case, the interchange forcing is relevant to the longer wavelengths k y p s < 
0.1, not to the range of largest linear growth rate. The reason is that the rms vorticity 
of the turbulence is strongly scale dependent, able to supersede the linear instability in 
the first decade longward of the p s scale, but quickly reduced in robustness at larger 
scale due to the factors of k± in the vorticity nonlinearity (note the factors of k± in 
Eq. 37). This is made clear if we plot the dynamical transfer spectra for the ExB vorticity 
against k± rather than k y , as in Fig. 17. The spectra are normalised to the total ExB 
energy. The linear growth stage (50 < t < 150) shows the linear instability peaking at 
k±p s ^0.S. In the transitional stage (200 < t < 250) the vorticity nonlinearity imposes 
its own character in that same spectral range. In the stage of fully developed turbulence, 
the larger scales find their saturated state, and in this case the interchange effects clearly 
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control the range k±p s < 0.1. Whether the turbulence is ultimately of the drift wave 
or interchange type depends on whether this larger scale range can decouple from the 
nonlinear vorticity dynamics which is mostly at 0.3 < k±p s < 1.0. If the linear effects 
are important only in this small scale range, then the linear instabilities are completely 
superseded by the turbulence and do not play a role in its dynamics. 

The linear analysis shows that ujb is often the dominant effect in creating a nonzero 
phase shift a p ^ and hence the linear instability, giving it a ballooning physical (not just 
morphological) character. But as the linear eigenmode reaches finite amplitude and the 
nonlinearities emerge, the turbulence not only saturates but changes its physical character. 
The mechanism maintaining the finite level of the parallel current changes from interchange 
forcing to the polarisation drift. We saw in the vorticity dynamical spectra in Fig. 6, that 
in the fully developed turbulence the interchange forcing does not assume control until the 
collisionality reaches what we found as the resistive ballooning threshold. The threshold 
is at about 

(^s)crit = (CW)crit = 1 (38) 

which is the boundary beyond which the interchange effects overcome the adiabatic re- 
sponse. At that threshold the transport does not sharply change, but the mode structure 
does so. So we can properly think of this transition as an eigenmode regime boundary, 
and of the turbulence to the low collisionality side of it as drift wave turbulence. 

The supersession of the linear instability spectrum found by the preceding analysis is 
the same as was recently shown in a detailed study of the drift wave nonlinear instability 
[39]. Not only does this form of self sustained turbulence cause turbulence in a sheared 
slab geometry when initialised at finite amplitude [10], it can also impose itself in favour 
of the linear instability mechanisms provided these are active in the same spectral range. 
The vorticity dynamics is active principally in the spectral range 0.3 < kj_p s < 1.0. Linear 
instabilities in this range are only relevant if they are strong enough to overcome the 
dynamical rate indicated by the rms vorticity of the turbulence — a reasonable rule of 
thumb would be a diamagnetic level of vorticity given by 

fl D = v D /L p (39) 

(unity in normalised units), where vd = c s p s /L p is the diamagnetic velocity. ExB flow 
shear must be stronger than this to suppress the turbulence [33] . Short wavelength linear 
instabilities will be relevant only if 

1l > (40) 
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i.e., their growth rate overcomes this. Long wavelength linear instabilities are always 
relevant because they can easily take over the spectrum against what is in this range a 
relatively weak vorticity nonlinearity. The result is, though resistive ballooning is not 
important to tokamak edge turbulence, ideal ballooning and other such large scale MHD 
instabilities are always relevant. 

The result that the turbulence changes character from the initial linear eigenmode due 
to supersession by vorticity scattering (O rms > 7^ but with little direct dissipation) is also 
interesting because the resulting nonlinear drive rate r + is quite smaller than the original 
7i, even for this borderline drift wave/interchange case. It shows that contrary to usual 
practice (cf. Ref. [23]), the judgement as to whether a linear or nonlinear instability has 
principal relevance should not rest on which one in its own native model produces the larger 
anomalous transport coefficient (i.e., net turbulent flux). If a simple comparison is to be 
made, the one between O rms and 7^ is more useful. But one should be careful about the 
/cj_-dependence of both, and ultimately a comprehensive set of numerical solutions to the 
nonlinear equations containing both models is required to make an informed judgement. 
The same is true of the nonlinear saturation mechanism, to which we now turn. 



VIII. Drive and Saturation Mechanisms 

The roles of the principal nonlinear ities, we • VVj_0 and • Vp e , are the same in 
toroidal geometry as in pure slab drift wave turbulence. We can show this by comparing 
the DALF3 model to the two dimensional drift wave model in a simple test in which we 
remove either of the two principal nonlinearities from the model. As a baseline case we 
take the DALF3 model run to t = 1000, well into saturation, for the nominal parameters 
given in Eq. (18), but with v = 2 (hence C = 5.1 and z/g = 0.25) and jUy = 0. Then, we 
go back to t = 500 and restart with either • VVj_</> removed ("novor") or all vg • V 
nonlinearities except v# ■ VVj_0 removed ("voronly") or vg- Vp e removed, and run again 
to t = 1000. The results are shown in Fig. 18. 

The vorticity nonlinearity very strongly excites the turbulence; it does not saturate 
even for Dirichlet boundary conditions, as without v# ■ Vp e there is no nonlinear mixing 
of the pressure. The pressure gradient continues to feed free energy into the disturbances, 
with the vorticity nonlinearity maintaining sufficient V|| Jy for the turbulence to access 
regions of statistical phase space with substantial phase shifts. The pressure nonlinearity 
on the other hand does not lead to arbitrarily strong turbulence, but saturates the overall 
amplitude at a somewhat higher level. The linear instability is back at work in the absence 
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of the nonlinear polarisation drift, balancing Vy Jy with K.{p e ). This result obtained for 
each of the five cases for which it was checked (u = 1, 2, 5 for the DW cases and v = 2, 5 for 
the BM cases); in no case does the turbulence saturate through the v#- VVj_0 nonlinearity. 

There are subtle differences to this result in the two dimensional (2D) systems. Though 
most of the subgrid dissipation in all the models goes through • Vp e , it remains true 
for the 2D Hasegawa-Wakatani model (where — C _1 V| is replaced by the constant D) 
that for many parameter combinations both nonlinearities are required in order to achieve 
saturation. Nevertheless, the strongest growth is obtained when v# • VV^_</> is present 
but \e ■ Vp e is absent. The main restriction of the 2D Hasegawa-Wakatani model is that 
V| cannot change its characteristic magnitude in response to changes in amplitude of the 
turbulence. 

The 2D model in which the nonlinear instability was first demonstrated is the the 2D 
sheared-slab model, a simple generalisation of the Hasegawa-Wakatani model to incorpo- 
rate magnetic shear: the collisional, electrostatic limit is obtained by taking (3 = fi = 0, 
neglecting any dependence in the s-direction, and replacing Vy with x(d/dy), that is, the 
coordinate system is not field aligned [10]. In that model the details of the dynamics are 
complicated by the way fey is tied to x; among the nonlinearities, retaining only ve • VVj_0 
causes the entire x-domain of the turbulence to be narrowed (as is the case with adiabatic 
electron turbulence [9]), while retaining only v# • Vp e causes the x-domain to spread. 
These results were found by taking the 2D slab restrictions within the present DALF3 
code and running with a resolution of 128 x 512 grid nodes in x and y, within the same 
xy-domain size of 20 x 80 in units of irp s , respectively. The cases run were C = 1, 10, and 
100. The effects on the morphology due to the absence of one of the nonlinearities was the 
same as in the 3D model: v# • VVj_0 produces isotropic, monopolar vortices in <p which 
excite p e unevenly but to large amplitude. In this case however, vg • Vp e by itself can 
produce the nonlinear instability at much the same level as with both nonlinearities, by 
acting indirectly through the partial time derivatives on both state variables. These details 
of nonlinear and drive and saturation deserve their own study (to be produced elsewhere); 
it suffices here to have presented the broad outlines. 

In three dimensions the degrees of freedom represented by x and fey are independent, 
and the situation is clearer: the slab model which one obtains by setting ujb = produces 
the same results [33] both quantitatively and qualitatively as in the toroidal case shown in 
this Section, except for the effect of the toroidal case's linear instability in the absence of 
Ve- VVj_0. Removal of the vorticity nonlinearity leads to linear drive (in the toroidal case) 
and saturation. Removal of the pressure nonlinearity leads to ever growing turbulence and 
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no saturation. The result is that in three dimensions the vorticity nonlinearity always 
drives and the pressure nonlinearity always saturates. One important consequence is that 
all scenarios which rely on Kelvin-Helmholtz effects on turbulent flows for saturation of 
tokamak edge turbulence are ruled out. 



Here we briefly confirm that the foregoing does not depend on the fact that we re- 
stricted to a model with isothermal electrons and cold ions. That was done to make the 
physics of the system as transparent as possible. To put the temperature dynamics back 
in, we refer to the methods of Ref. [2] to treat the parallel dynamics; the heat fluxes are 
themselves dynamical variables so that the two pairs of additional dependent variables 
{T e ,5e||} and {Ti,qi\\} are treated on the same self consistent footing as {n e , Jy} in the 
DALF3 model (in which p e = n e ). This treats time dependent Landau damping, extending 
the model to the weakly collisional regime in which the Braginskii equations lose validity 
(mainly by overestimating damping of T e and especially Ti, since the turbulence is about 
two orders of magnitude faster than the ion collision frequency). We now have p e = n e + T e 
and pi = Ti?i e + Ti as normalised pressure disturbances, where Ti = Ti/T e gives the warm- 
ness of the ions. The only remaining subtlety is gyroviscosity. The form of the diamagnetic 
cancellation used here is 



in physical units, where u* is the ion diamagnetic velocity, 11* is the gyroviscosity tensor 
(the diamagnetic momentum flux), and x is a scalar involving finite gyroradius effects 
and divergences of both the heat flux and the ion velocity [40]. The sole effect of the 
gyroviscosity is to cancel u* in the advection. The resulting normalised equations are [41] 



IX. The DALFTI Model and Results for Warm Ion Turbulence 
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dun _ _ 2 _ 

= - V|| {Pe +p t +Pe+ Pi) + M|| V||U|| (47) 

A^ + aLe(5e||) = -|V||(T e +f e ) - ^/^ e (§e|| +« e J||) (48) 

e-^- + a Li (qq) = -^Vy (T, + T;) - ^ez/^ y (49) 
with current, pressures, and temperature ratio 

jjl = -Vx^4|| p e = n e + f e pi = nn e +fi n= Ti/T e (50) 

and Landau damping operators 

a Le = jj 1/2 (l - 0.2Vjj) a Ll = (ne) 1 ' 2 (l - 0.2V|) (51) 

including linear gradient terms through 

n e = -u n x T e = -oj t x Ti = -nuix p e = -(oj n + u t )x pi = -Ti{u n + Ui)x 

(52) 

The geometry is the same as described in Section II. This is called the DALFTI model, 
extending drift Alfven turbulence to the warm ion regime with Landau electrons and ions. 

The extra parameters describing collisional dissipation are the normalised collisional 
frequencies, 

= Vi=^~ (53) 

where the r ve ^ are the standard Braginskii collision times [30], in which terms the drift 
wave collisionality is C = r\pi,v. The numerical constants are 

rj = 0.51 a e = 0.71 K e = 3.2 Ki = 3.9 (54) 

for pure hydrogen, treating resistivity, the thermoelectric effect, and electron and ion 
thermal conduction, respectively. With these numbers C has the same meaning as in the 
DALF3 model. If the v e ^ become arbitrarily large we recover the Braginskii equations, as 
discussed in Ref. [2], but we are not in that regime for typical tokamak edge parameters, 
especially for ions. 

We briefly report the reconsideration of the cases in Section V for the DALFTI model, 
with the same nominal parameters as in Eqs. (18), for which v = 3, and with u n = ut = 
Ui = Ti = 1 placing the model in the drift wave gradient regime (ion temperature gradient 
turbulence, with larger values of uji/u n , will be treated elsewhere; cf. also Refs. [33] and 
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[41]). The warm ion MHD model is formed analogously to DALF3: the combinations 4>+Pi 
and (f) — p e are all replaced simply by 4> (including the background gradient pieces where 
they appear under w E ■ V and Bj_ -V). In Fig. 19 we find the transport scaling and the 
relative disturbance amplitudes for cases with various z/, with the four main drift wave cases 
with v = 1, 2, 5, 10 checked with doubled resolution in the drift plane. These results may 
be compared to Fig. 3. Two things should be noted: The results are somewhat clouded 
by the presence of the ion temperature disturbance, T^, which always exhibits interchange 
dynamics due to the fact that the nominal sound wave transit frequency c s /qR is typically 
smaller than the dynamical frequencies in the range of 0.2c s /Lj_, so while the interchange 
effects are relatively weak compared to the turbulence they are still much stronger than the 
parallel ion dynamics. We also have the fact that Ti plays no role in the parallel electron 
dynamics, and is the only state variable for which this is true. Secondly, p e includes both 
n e and T e , so that the ratio <j>/n e is also influenced by changes in the relative role of 
the temperature and density gradients. This follows from the fact that more nonadiabatic 
electron dynamics in the stronger turbulence for higher collisionality leads to comparatively 
stronger effect due to VT e [2,10]. Nevertheless, although the transition regime is wider 
than in cases with pure electron dynamics (cold ions), it still occurs for vb only slightly 
less than unity. 

The transition to resistive ballooning for these warm ion cases is shown in Figs. 20- 
22. As v is swept through the values 2, 5, and 10, we find the gradual increase in the 
amplitude of 4> relative to the other state variables. In the drift wave regime T e roughly 
has the same structure as h e in the DALF3 model (compare with Fig. 9); indeed, just as in 
slab drift wave turbulence these two quantities control the release of free energy from the 
electron thermal gradient and hence play the same role in the turbulence [10]. The fact 
that they are less ballooned than the density or ion temperature (or the transport fluxes) 
is the reason that drift wave mode structure persists in turbulence in toroidal geometry, 
just as for h e in the DALF3 model. The arrival of the pure MHD regime for v = 10 is 
also facilitated by the contribution of pi to the total ion flow stream function, W — (f> + Pi, 
and hence the vorticity, V 2 ± W. This adds to the power of the total nonlinear polarisation 
drift and removes the linear forcing terms from dominance at lower collisionality than in 
a pure MHD model for which the ion flow is just the ExB flow (cf. also [1]). Tests have 
shown that the presence of the nonlinear gyroviscosity has a large effect in the practical 
drift wave regime 1 < v < 5, but we leave the details of this to a subsequent publication 
(without the gyroviscosity, the vorticity in Eq. 42 is advected by the total flow expressed 
with W, and there is no correction; in contrast to the actual situation, both energy and 
enstrophy are conserved in the case that gyroviscosity is neglected). Due mostly to the 



26 



new effects brought in by the finite ion temperature, the onset of the transition to the 
resistive ballooning regime is moved to about 

Merit = 0.5 (55) 

that is, downward by about a factor of two, generally reflecting the equal and additive 
effects of the electron and ion pressures. 

X. The Ideal Ballooning Boundary 

Up to now we have focussed on the comparison between drift wave and resistive 
ballooning turbulence, where only the collisionality is varied. The other boundary is better 
known: the ideal MHD boundary, which gives the onset of dissipation free ballooning 
modes called ideal modes [15]. Experimentally this is thought to give the beta limit, which 
is known both globally and locally [42] . In Fig. 23 we show this for both the DALF3 and 
DALFTI models, setting v = 1 to get a sharper transition. The transport is computed as 
a flux as before, but this time it is given in terms of the traditional transport coefficients 
in physical units - - ambiguity with regard to the trend is minimised by choosing all 
gradient scale lengths equal to L± and by avoiding the effect of hidden normalisation, as 
the diffusivities scale with p 2 s c s /L±_ and hence with (/3/C) 1 / 2 if B and the scale lengths 
are all held fixed. We find a clear boundary in both models though it is accentuated by 
the presence of T^. The jump in the transport appears at different values of (3 because the 
total pressure differs in the models. The standard ideal MHD stability parameter is 

«m = -q 2 RV(3 = (3uj B [(uj n + u> t ) + Ti(uj n + u>i)} (56) 

in both physical and normalised units, where the two combinations in parentheses give 
the normalised Vp e and Vpi, respectively, giving a factor of 4 for these parameters for 
DALFTI and a factor of 1 for DALF3. The transport amplitude is plotted against ctM, 
showing that the effective ideal ballooning boundary for the turbulence is given by 

0.2 < (a M )crit < 0.6 

for s = 1. The regime change starts closer to 0.2 but the mode structure changes are 
not complete until the longer wavelengths separate from the rest of the dynamics, for «m 
about 0.6. 

The mode structure changes in the ideal ballooning transition are shown in Figs. 24- 
26, for the three cases f3 = {0.3, 1.0,3.0} for the DALFTI model, measured in the same 
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way as in Section V. The signature of the change in character of the turbulence is very- 
clear in these figures. When the ideal ballooning takes over, it does so in the longest 
wavelengths as their phase shifts separate from the rest of the spectrum. At the same 
parameters, the cross coherence between n e and <fi is lost. The signatures of the change 
in character of the turbulence in this transition is sharper and clearer than in the one for 
resistive ballooning. All vestige of drift wave character in the cross coherence and parallel 
structure diagnostics is lost when the ideal ballooning boundary is crossed. The phase 
shift diagnostic shows that the long wavelength piece splits off and is no longer influenced 
by the vorticity dynamics coming from the shorter wavelengths, with the result that the 
entire spectrum no longer acts as a self consistently causal unit. The same behaviour is 
exhibited by both the DALF3 and DALFTI models. 

The ideal ballooning cases saturate only by depleting the pressure profile: With pe- 
riodic boundary conditions in x these cases do not saturate at all, instead forming a 
wide jet in the x-direction with k x = and k y L x rj n where L x = 2ir/AK is the x- 
domain width. With Dirichlet boundaries the jet becomes a more or less isotropic cell 
with k x rs k y as tt/L x , saturating only as the resulting flux balances the source terms 
maintaining the profile in the edge damping zones. This means that saturation proceeds 
only directly on the profile, with the combined action of the nonlinearities unable to suf- 
fice. In no case was the vorticity nonlinearity observed to provide saturation through 
Kelvin-Helmholtz effects. This has been found to be a feature of the "thin atmosphere" 
situation with aspect ratio A >> 1, contrary to other efforts [4] which use A = 1, artificially 
constricting the formation of wide down-gradient flows. 

XI. Summary — Turbulence in Context 

What these results show is that it is important to consider the turbulence in context 
when making judgments about its character, especially when in a situation like this there 
is more than one possible mode structure into which it can arrange itself. Moreover, it 
is essential to set up computations with equal regard for the properties of both or all 
the possible eigenmode types, so that the results are not prematurely anticipated. We 
find that this physical character undergoes strong changes no only in spatial structure 
but also in energetics as the linear instability makes the transition into turbulence. Since 
the turbulence has no detailed memory beyond several (rather short) correlation times, 
the character the dynamics has in the linear stage is not only not relevant, it delivers an 
incorrect paradigm as to understanding of the basic nature of the turbulence. This serves 
to underscore the danger of relying on linear instability theory in the formative stages of a 
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body of work whose aim is to understand turbulence and transport. The linear instability 
is relevant whenever it can act on the longer wavelengths (low k±), where the vorticity 
nonlinearity is weak, but at high k±p s ~ 1 (i.e., "high-n ballooning") they lose relevance 
because these short wavelength instabilities are superseded by the self sustained drift wave 
turbulence. Since this turbulence has an rms vorticity greater than the ideal interchange 
growth rate, it scatters the linear eigenmodes apart before they can grow (in analogous 
fashion to the way a background ExB vorticity suppresses instabilities) . 

The important thing to note about these results is that they all speak together; there 
are no contradictions. In all cases, the mode structure for low collisionality and low 
beta exhibits clear drift wave character, and the transition to resistive MHD is found 
to start for v > 3, which for this choice of parameters corresponds to = 41, where 
= (gi?/r iye V e )e _3 / 2 and in normalised units = z/ / u 1 / 2 e -3 / 2 or about 13.6z/ e and 5.3C 
for the nominal scale ratios in Eq. (18), having assumed an inverse aspect ratio of e = 0.3. 
Therefore, having examined the properties of the turbulence in its native context, we 
can determine where the MHD model begins to be valid, and where we need the drift 
wave model not only for computation but for fundamental understanding. Where the 
MHD model loses validity due to the adiabatic response becoming important, it also loses 
usefulness as a paradigm, because the basic physics of how the disturbances in pressure 
and ExB flow communicate undergoes fundamental changes. Generally, MHD character 
is found when the interchange forcing on the pressure overcomes the adiabatic response. 
For ideal MHD the regime boundary is the ballooning limit at the critical value of olm = 
q 2 R |V/3| for the particular geometry; in the circular tokamak model it is near (ctM)crit = s. 
For resistive ballooning the ideal Alfven response is replaced by its dissipative MHD limit, 
and the boundary is near vb = Cub = 1- The factors of order unity will vary with magnetic 
geometry and some of the complications of the temperatures, but generally for the low beta 
and moderate collisionality regime in which fusion plasmas are found, one either finds drift 
wave turbulence if the equilibrium is steady, or ideal ballooning phenomena when there 
are disruptive events which typically collapse the pressure gradient. 

Concerning the situation of linear instabilities vis-a-vis self sustained drift wave tur- 
bulence, we find a general indication that short wavelength (k y p 8 ^0.3) linear instabilities 
are relevant to the turbulence if and only if their linear growth rates are stronger than 
the general diamagnetic vorticity level: 7^ > Qd. Long wavelength (k y p s ^ 0.1) linear 
instabilities, on the other hand, are always relevant because the vorticity nonlinearity is 
relatively weak in that spectral range. Borderline cases can only be decided with well 
resolved computations assisted by detailed diagnosis. 
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It is important to note that the question of drift wave or ballooning character is one 
which is decided by the physical processes responsible for free energy generation, transfer, 
and saturation, not by the general look and feel of the turbulence which would be a 
very arbitrary issue. The formative literature for both eigenmode types made very clear 
statements of what these processes should be, and in order to decide which is most relevant 
one must diagnose those processes directly and within the context of the turbulence itself. 
It is precisely that which we have done herein. 
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Figures 




Figure 1. (left) Time traces of the half squared amplitudes of p e and 4>, re- 
spectively labelled by A v and A^, in the drift wave model (right) and the 
interchange model (right), showing saturation. Only in the drift wave model 
are the time traces for A p and A^ similar. 
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Figure 2. (left to right) Drive and transfer spectra, and cross coherence and 
phase shifts between p e and 0, for the drift wave (top row) and interchange 
(bottom row) models. In the drift wave case the sink spectrum (labelled 'C') 
is relatively flat, and the transfer ('tr') is also due to V\\J\\ and is positive at 
short wavelength and negative at long wavelength, while in the interchange 
case it is due to K,{p e ) which follows the source spectrum ('pe'). The drift wave 
case shows strong cross coherence and a narrow phase shift distribution closer 
to zero than to tt/2. The interchange case shows dominance by the longest 
wavelengths, no cross coherence, and phase shifts near tt/2 due to the strong 
driving and weak coupling. The phase shift distributions contours are 0.3, 0.5, 
and 0.8 times the maximum, and for the cross coherence they are 0.37, 0.5, and 
0.8 times the maximum. 
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Figure 3. (left) Transport scaling of drift wave and interchange turbulence in 
toroidal geometry, from the DALF3 and reduced resistive MHD models labelled 
'DW and 'BM', respectively. The drift wave cases show the clearest scaling with 
collisionality at low C = 2.55z/. At asymptotically large C the trends will merge, 
but that limit is not reached, (right) Half squared amplitude ratio (including 
only k y 0) for the two sets of cases. Due to the adiabatic response, <p tracks 
p e for drift wave turbulence, but in the MHD model of the interchange cases p e 
is unaffected by the Alfven dynamics, and so instead of 4> being forced towards 
p e it is forced towards zero. The extra points marked with squares for DW 
{v = 1, 2, 5, 10) are with double resolution in the drift plane. 
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Figure 4. Cross coherence between p e and <p, for drift wave (top row) and in- 
terchange (bottom row) turbulence, for v = 5, 10, and 20 (left to right), where 
C = 2.55u and ug = Cub as defined in Eq. (32). Compare with the results 
in Section IV for drift wave and interchange turbulence. The DALF3 model 
results in drift wave mode structure even for larger z/g, while the MHD model 
always shows interchange mode structure. Contours are as in Fig. 2. 
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Figure 5. Phase shift distributions of p e ahead of <f) at each k y , for drift wave 
(top row) and interchange (bottom row) turbulence, for v = 5, 10, and 20 (left 
to right), where C = 2.55z/ and v& = Cub as defined in Eq. (32). Compare 
with the results in Section IV for drift wave and interchange turbulence. The 
DALF3 model results in drift wave mode structure for v < 10, while the MHD 
model always shows interchange mode structure. The transition to resistive 
ballooning turbulence in the DALF3 model occurs in the longest wavelengths, 
kyp s < 0.1. Contours are as in Fig. 2. 
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Figure 6. Spectra of the rms transfer dynamical levels for each k y in the spec- 
trum, comparing the sizes of 0V|| J\\ ('j'), (j)lC(p e ) ('k'), and 4>ve • W\4> ('e'), 
for drift wave (top row) and interchange (bottom row) turbulence, for v = 5, 
10, and 20 (left to right), where C = 2.55z/ and z/g = Cub- In drift wave tur- 
bulence the transfer through the current is too large to be accounted for by the 
curvature, and at all wavelengths the vorticity nonlinearity is balanced only 
by the linear time derivative. Nonlinear vorticity dynamics is generally much 
stronger when the adiabatic coupling mechanism p e <-> Jy is present, leading 
to a self consistent situation in which both mechanisms catalyse each other. 
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Figure 7. Amplitude spectra of p e , 0,and V^0, respectively labelled by 'pe', 
'phi', and 'vor', for drift wave (top row) and interchange (bottom row) turbu- 
lence, for v = 5, 10, and 20 (left to right), where C = 2.55z/ and v B = Cujb- 
With only the density present one cannot distinguish the mode structure or 
dynamics, but if <fi is also present the amplitude ratio, shown in Fig. 3, is de- 
cisive, as the spectrum of <p follows that of p e only for drift wave turbulence. 
The transition to resistive ballooning turbulence in the DALF3 model {<p ^> p e ) 
occurs in the longest wavelengths, k y p s < 0.1. 
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Figure 8. Spectra of the ExB gradient drive ('pe'), dissipation ('C'), and the 
magnetic flutter drive ('mag'), for drift wave (top row) and interchange (bottom 
row) turbulence, for v = 5, 10, and 20 (left to right), where C = 2.55z/ and 
vb = Cub- For drift wave turbulence the shorter wavelengths contribute more 
to the energetic drive and hence the ExB transport. Compare the positions of 
the energy containing range (Fig. 7), the energy producing range (this figure) 
and the vorticity catalysing range (Fig. 6). These show that for drift wave 
turbulence the entire spectrum acts as a single unit. 
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Figure 9. Mean squared amplitude envelopes (k y 7^ only) showing parallel 
structure of cf>, p e , and h e = p e — (f), respectively labelled by 'phi', 'pe', and 'he', 
for drift wave (top row) and interchange (bottom row) turbulence, for v = 5, 
10, and 20 (left to right), where C = 2.55z/ and v B = Cujb- Drift wave mode 
structure is exemplary for v < 10 in the DALF3 model, with h e flatter and 
smaller than either p e or 0. 
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Figure 10. Amplitude envelope of the ExB transport ('pe'), for drift wave (top 
row) and interchange (bottom row) turbulence, for v = 5, 10, and 20 (left to 
right), where C = 2.55z/ and v B = Cujb- The magnetic flutter transport is 
negligible on this scale. The ballooning in the transport becomes somewhat 
more pronounced in the transition to interchange character for v > 10, but for 
all cases the ballooning is much stronger for the MHD model than for DALF3. 
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Figure 11. Time evolution of the DALF3 case with v = 10 (hence C = 25.5 and 
vb = 1.25) out of a linear initial state at small amplitude, run to t = 1582. (top 
left) Half squared amplitudes of <j) (denoted A^) and p e (A p ). The larger ampli- 
tude departures of A^ from A p reflect zonal flow activity in the fully developed 
turbulence, (top right) The transport caused by the turbulence, (bottom left) 
Growth rate (IV), gradient drive rate (r+), and total dissipation rate (Fe)- 
(bottom right) Vorticity (rms) compared to r_|_ and Ft- The linear mode leads 
to overshoot, initial saturation (t ~ 200) is reflected in the first drop of Ft 
to zero, and then the nonlinear mode structure takes over shortly thereafter, 
with full development with robust transport and zonal flow activity established 
after t = 400. The turbulence not only saturates but changes character when 
the rms vorticity overcomes the linear growth rate ('Jl, equal to Ft in the lin- 
ear stage). Nonlinear saturation of the linear instability would obtain if n r ms 
were comparable to 71,, but the situation with O rms 7^ indicates complete 
supersession of the instability by turbulence which has its own dynamics. 
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Figure 12. Saturation and initial transition to turbulence, as seen in the spatial 
morphology of <fi an d p e in the linear growth stage to t = 190, and the initial 
saturation stage after t = 200. The transition between linear and nonlinear 
mode structure is most clear in the disappearance of x-direction flows rem- 
niscent of bouyant plumes. This represents a transition away from interchange 
dominated dynamics to a more isotropic turbulence as the vorticity nonlinearity 
replaces the interchange forcing as the principal mechanism supporting finite 
parallel currents (nonadiabatic electron dynamics). (Positive/negative values 
are indicated by solid/dashed lines. Only 1/4 of the computational domain in 
the y-direction is shown. The exact moment of saturation is t = 201, and the 
growth rate becomes positive again at t = 212.) 
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Figure 13. Transition to fully developed turbulence, as seen in the spatial mor- 
phology of 4> and p e through the stage of nonlinear structure adjustment, 
and then in the stage of statistical saturation in which the nonlinear growth 
rate fluctuates near zero. As the saturated state finds itself, the scale of mo- 
tion increases, and until t = 345 the nonlinear growth rate is positive. Zonal 
flows emerge after about t = 400 and reach statistical equilibrium after about 
t = 600. The zonal flows are part of the nonlinear mode structure, but by 
the time they emerge the interchange driven flows of the linear stage are long 
gone. The correlation time is about 6 in these units. (Positive/negative values 
are indicated by solid/dashed lines. Only 1/4 of the computational domain in 
the y-direction is shown. The exact moment of saturation is t = 201, and the 
growth rate becomes positive again at t = 212, and negative again at t = 345.) 
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Figure 14. Saturation and transition to turbulence, as seen in the dynamical 
transfer spectra for the ExB vorticity, for the linear stage to about t = 190, 
the saturation stage around t = 200, and the nonlinear structure adjustment 
stage after about t = 210. Compare with Fig. 6. Initially the eigenmode is 
controlled by interchange forcing, with the linear polarisation drift negligible 
(cf. Eq. 35). But as the amplitude becomes finite the vorticity nonlinearity, the 
same one which causes the drift wave nonlinear instability, emerges to become 
the principal agent supporting nonadiabatic electron dynamics. The transition 
is extremely rapid, taking place within about one correlation time for the fully 
developed turbulence. (The exact moment of saturation is t = 201, and the 
growth rate becomes positive again at t = 212, and negative again at t = 345.) 



46 





Figure 15. Saturation and transition to turbulence, as seen in the parallel enve- 
lope structure (k y ^ only), in the linear stage to about t = 190, the saturation 
stage around t = 200, and the nonlinear structure adjustment stage after about 
t = 210. Compare with Fig. 9. The principal signature of the transition is the 
emergence of <p supported by the vorticity nonlinearity. The interchange flows 
of the linear stage are eliminated by the turbulent vorticity, and the adiabatic 
response causes <j> to track p e . The degree of asymmetry is also reduced, es- 
pecially in h e . This transition between linear (h e ~ p e ^> 0) and nonlinear 
(4> ~ p e 3> h e ) mode structure is almost as rapid as in Fig. 14. (The exact 
moment of saturation is t = 201, and the growth rate becomes positive again 
at t = 212, and negative again at t = 345.) 
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Figure 16. Phase shift distributions of p e ahead of <fi at each k y , for the lin- 
ear stage, averaged over 50 < t < 150, through saturation, averaged over 
200 < t < 244, and for the stage of fully developed turbulence, averaged over 
502 < t < 611. Compare with Fig. 5. The turbulence emerges to supersede 
the linear structure with its own, due to the fact that the rms vorticity level 
of the turbulence is larger than the linear growth rate of the original instabil- 
ity. The linear instabilities in the range 0.3 < k y p s < 1.0 have no role in the 
turbulence. Only in the fully developed stage does the interchange dynamics 
for kyp s < 0.1 emerge to make this case with vb = 1-25 the transitional one 
between turbulence of the drift wave and interchange type. 
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Figure 17. Dynamical transfer spectra for the ExB vorticity, for the linear stage, 
averaged over 50 < t < 150, through saturation, averaged over 200 < t < 244, 
and for the stage of fully developed turbulence, averaged over 502 < t < 611, 
plotted against k± rather than k y , showing the scale of motion rather than 
the wavelength in the drift direction. The linear instability is dominantly in 
the same range which is later dominated by the turbulence, making the linear 
interchange dominated mode irrelevant. When interchange effects do enter, as 
in this transitional case with v B = 1-25, they do so at larger scale where they 
more easily overcome the native vorticity of the drift wave turbulence. 
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Figure 18. Saturation mechanism in drift wave turbulence in toroidal geome- 
try. The run with v = 2 (hence C = 5.1 and v B = 0.25) taken to t = 1000 
(solid curves) is restarted from t = 500 (dashed curves) with the ExB vorticity 
nonlinearity either left out ('novor') or with all ExB nonlinearities except the 
vorticity one left out ('voronly'). Without the vorticity nonlinearity, the lin- 
ear drive is balanced by mixing via vg • Vp e and saturation occurs. Without 
the pressure nonlinearity, the vorticity is vigorously scattered via w E ■ VVj_0 
with little dissipative effect, and this nonlinear excitation continues indefinitely 
without saturation. 
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Figure 19. Transport scaling (left two frames) and amplitude ratios (right) of 
drift wave (solid lines) and resistive MHD (dashed lines) turbulence in toroidal 
geometry under the DALFTI model. Both ion and electron heat transport are 
shown as a pair; in each case the ion transport curve is the one lying slightly 
higher in the pair. The MHD model, which neglects the drift wave coupling 
terms between p e and Jy, is insensitive to collisionality (C = 2.55u) because 
although <p is too small, pi is too large, compared to the drift wave model. The 
transition to the MHD regime at a lower v > 3 is assisted by the warm ion 
interchange physics, since Tj does not feel the adiabatic response. The extra 
points marked with squares for DW (u = 1, 2, 5, 10) are with double resolution 
in the drift plane. Compare with Fig. 3. 
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Figure 20. Cross coherence between n e and 0, for drift wave (top row) and 
interchange (bottom row) turbulence, for v = 2, 5, and 10 (left to right), where 
C = 2.55z/ and vb = Cub as defined in Eq. (32). The turbulence makes the 
transition from drift wave to resistive ballooning mode structure for vb of about 
0.5 (here, v = 5). Compare to Fig. 4. Interchange turbulence with warm ions 
is very violent at small scales due to the effects of gyroviscosity, reducing the 
timestep and shortening the runs. 
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Figure 21. Phase shift distributions of n e ahead of <f> at each k y , for drift wave 
(top row) and interchange (bottom row) turbulence, for v = 2, 5, and 10 (left to 
right), where C = 2.55z/ and v B = Cujb as defined in Eq. (32). The turbulence 
makes the transition from drift wave to resistive ballooning mode structure for 
vb of about 0.5. Compare to Fig. 5. The tendency of the phase shift to go 
to — 7T at small scales is the signature of the nonlinear ity in the gyroviscosity 
which conserves energy but not vorticity for warm ions, producing the effects 
seen in Fig. 20. As in the DALF3 model, the transition to resistive ballooning 
turbulence in the DALF3 model occurs in the longest wavelengths, k y p s < 0.1. 
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Figure 22. Mean squared amplitude envelopes (k y 7^ only) showing parallel 
structure of 0, n e , T e , and Tj, respectively labelled by 'phi', 'ne', 'Te', and 
'Ti', for drift wave (top row) and interchange (bottom row) turbulence, for 
v = 2, 5, and 10 (left to right), where C = 2.55z/ and v B = Cujb as defined 
in Eq. (32). The turbulence makes the transition from drift wave to resistive 
ballooning mode structure for vb of about 0.5, with additional effects due to the 
contribution of pi to the vorticity, as described in the text. Note p e = n e + T e 
and pi = Ti?i e + Ti in this model. 
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Figure 23. Transport scaling as a function of the ideal ballooning parameter ccm, 
for the DALF3 ('3') and DALFTI (T) models, expressed as diffusivities, with 
Xe,i including the convective contributions. The collisionality was C = 2.55, 
in the drift wave regime. The transport is given in physical units (m 2 /sec, 
assuming deuterium ions, and B = 2.5 T, and L± = 4.2 cm), compensating for 
the effect of varying (3 on the normalisation scale p^c s /L±. The ion temperature 
assists the transition to ideal MHD in the measure and for the same reason as 
for the resistive ballooning cases: more total pressure gradient, and T; > T e . But 
the ideal ballooning threshold, now between 0.2 and 0.6, is lowered by a factor 
of at least two relative to the linear analysis. The experimentally interesting 
range is actually «m > 0.2, corresponding to $ > 1 for the DALFTI model. 
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Figure 24. Cross coherence between n e and <f>, for the DALFTI model at C = 
2.55 hence vb = 0.25, for various «m in the range of the sharp transport rise, 
noting «m = 0.2(3. The turbulence makes the transition from drift wave to 
ideal ballooning mode structure in this parameter range. Compare to Fig. 20. 
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Figure 25. Phase shift distributions of n e ahead of <p at each k y , for the DALFTI 
model at C = 2.55 hence vb = 0.25, for various «m in the range of the sharp 
transport rise, noting = 0.2/3. The turbulence makes the transition from 
drift wave to ideal ballooning mode structure in this parameter range. Compare 
to Fig. 21. Here as well, the transition to resistive ballooning turbulence in the 
DALF3 model occurs in the longest wavelengths, k y p s < 0.1. 



57 




Figure 26. Mean squared amplitude envelopes (k y 7^ only) showing parallel 
structure of 0, n e , T e , and Tj, respectively labelled by 'phi', 'ne', 'Te', and 'TV, 
for the DALFTI model at C = 2.55 hence vb = 0.25, for various «m in the 
range of the sharp transport rise, noting olm = 0.2/3. The turbulence makes the 
transition from drift wave to ideal ballooning mode structure in this parameter 
range. Compare to Fig. 22. Note p e =n e + T e and pi = T(h e J rT i in this model. 
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